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ABSTRACT 

We report simultaneous millimetre and infrared observations of a sequence of 
very large amplitude quasi-periodic oscillations from the black hole X-ray binary GRS 
1915+105. These oscillations are near the end of a sequence of over 700 repeated events 
as observed at 15 GHz, and are simultaneous at the mm and infrared wavelengths to 
within our time resolution (< 4 min), consistent with the respective emitting regions 
being physically close near the base of the outflow. One infrared event appears to have 
no mm counterpart, perhaps due to highly variable absorption. The overall radio-mm- 
infrared spectrum around the time of the observations does suggest some absorption 
at lower frequencies. We calculate the energy and mass-flow into the outflow for a 
number of different assumptions, and find that the time-averaged power required to 
produce the observed synchrotron emission cannot be much less than 3 x 10 38 erg s -1 , 
and is likely to be much larger. This minimum power requirement is found regardless 
of whether the observed emission arises in discrete ejections or in an internal shock in 
a quasi-continuous flow. Depending on the similarity of the physical conditions in the 
two types of ejection, GRS 1915+105 may be supplying more power (and mass, if both 
have the same baryonic component) to the jet during periods of repeated oscillations 
than during the more obvious larger events. 
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1 INTRODUCTION 

The 'micro-quasar' GRS 1915+105 is one of the most cel- 
ebrated and widely-studied astrophysical objects of recent 
years. The system is extremely luminous and variable in 
both hard and soft X-rays (e.g. Foster et al. 1996; Morgan, 
Remillard & Greiner 1997; Belloni et al. 2000) and is a source 
of relativistic jets observed on arcsec and milliarcsec angular 
scales (Mirabel & Rodriguez 1994, hereafter MR94; Fender 
et al. 1999, hereafter F99; Rodriguez & Mirabel 1999, here- 
after RM99; Dhawan, Mirabel & Rodriguez 2000). Sams, 
Eckart & Sunyaev (1996) have reported extended infrared 
emission from GRS 1915+105, but its relation to the radio 
ejections is at present unclear. 

X-ray dips on timescales of minutes have been inter- 
preted by Belloni et al. ( 1997a, b) as the repeated disap- 
pearance and refill of the inner accretion disc, possibly due 
to extremely rapid transitions between 'canonical' black hole 
accretion states (Belloni et al. 2000). Pooley & Fender (1997; 
hereafter PF97) reported radio oscillations associated with 



such dips, and Fender et al. (1997; hereafter F97) discovered 
infrared analogs of these oscillations. The flat spectrum and 
correlated radio : infrared behaviour suggested that non- 
thermal synchrotron emission extended from the radio to 
the infrared regimes, the first time such high-frequency syn- 
chrotron emission had been observed from an X-ray binary 
(F97). Combined with the unstable accretion disc model 
of Belloni et al. (1997a,b) we suggested that a fraction of 
the inner disc was being repeatedly accelerated and ejected 
from the system (F97; PF97). Eikenberry et al. (1998a) con- 
firmed the association between X-ray and infrared events, 
and Mirabel et al. (1998, hereafter M98) clearly observed the 
correlation between X-ray, infrared and radio behaviour in 
the source. Additional simultaneous observations (Fender & 
Pooley 1998, hereafter FP98) showed a very clear correlation 
between sequences of oscillations at radio and infrared wave- 
lengths. Delays between different radio bands (PF97; M98) 
and between the infrared and radio bands (M98; FP98) 
clearly indicate that optical depth effects play an important 
role in the observed emission from these ejections. 
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Figure 1. Simultaneous mm (1.3 cm) and infrared K-band (2.2 fim) light curves of GRS 1915+105. The infrared data have been 
dereddened by Ak = 3.3 mag. These are the largest repeated oscillations ever observed from GRS 1915+105. 



Eikenberry et al. (1998b) showed that infrared emission 
line strengths vary in proportion to the continuum during 
oscillations. A single X-ray dip, spectrally associated with 
the temporary disappearance of the inner accretion disc, was 
also found to coincide with a small radio flare (Feroci et al. 
1999). More recently Eikenberry et al. (2000) report faint 
infrared flares whose association with the X-ray behaviour 
is uncertain, and Ogley et al. (2000) report significant flux 
from GRS 1915+105 at sub-millimetre wavelengths. 



2 OBSERVATIONS 

GRS 1915+105 was observed simultaneously on 1999 May 20 
with the United Kingdom Infrared Telescope (UKIRT) and 
the James Clerk Maxwell Telescope (JCMT), both located 
on Mauna Kea, Hawaii. 



2.1 UKIRT 

GRS 1915+105 was observed with IRCAM3 in the infrared 
K-band (2.2/im) on 1999 May 20, simultaneously with the 
longer duration of JCMT SCUBA observations (see be- 
low). Data reduction and calibration were performed with 
iraf, along the lines described in F97. Five clear oscilla- 
tion events were detected. The undereddened infrared flux 
densities reached 25 mjy at the peak of the oscillations, 
the largest amplitude oscillations reported to date in the in- 
frared. The data are plotted in Figs 1 & 2, dereddened by 
Ak = 3.3 mag (F97; this value is still rather uncertain). 
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Figure 2. Enlargement of the simultaneous mm/infrarcd obser- 
vations. While four of the five infrared oscillations correlate well 
with the mm emission, one event is clearly anomalous, with little 
mm response to the second infrared oscillation. 



2.2 JCMT 

The 1350 jj,m detector of the SCUBA instrument (Holland 
et al. 1999) on JCMT was used in the photometry. Each 
integration lasted approximately 4 min. Calibration of the 
flux-density scale used observations of Mars and Uranus. 
The airmass ranged from 1.01 to 1.92 during the observa- 
tions, and the optical depth at 1350 /im was less than 0.2 
throughout. The data are plotted in Figs 1 & 2. 
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2.3 Radio 

In order to piece together the composite radio-mm-infrared 
spectrum of the source at the epoch of our observations, 
we have utilised radio data from two different monitoring 
programs. Firstly we have used public data at 2.3 & 8.3 
GHz from the Green Bank Interferometer (GBI) monitor- 
ing program (e.g. Waltman et al. 1994). Observations at 15 
GHz with the Ryle Telscope (RT, e.g. PF97) reveal strong 
oscillations for at least 9 days before and 2 days after the 
simultaneous UKIRT/JCMT observations (Fig 3), with a 
slowly rising trend in mean level and amplitude. 



3 TEMPORAL AND SPECTRAL BEHAVIOUR 

The JCMT observations reveal a sequence of 15 millimetre- 
wavelength oscillations with a quasi-period of ~ 23 min. 
Four of the oscillations have been observed simultaneously 
in the near-infrared, and both the (dereddened) near in- 
frared and mm oscillations have an amplitude of 300-350 
mjy. These are by far the largest amplitude oscillations ever 
observed from GRS 1915+105, including radio wavelengths 
(see e.g. PF97 for 'typical' radio oscillations) . There is no de- 
tectable delay between the emission at the two wavelengths, 
to the time resolution of the JCMT data, ~ 4 min. 

F97, M98 and FP98 all report infrared oscillations of 
comparable amplitude (when dereddened by Ak =3.3 mag) 
to radio oscillations observed around the same time. Al- 
though radio observations were not made strictly simultane- 
ously with these mm/infrared observations, monitoring with 
the RT clearly reveals radio oscillations, with a comparable 
quasi-period of ~ 20 min, to have been occurring for at least 
two days before and afterwards (Fig 3). However, the ampli- 
tude of the mm and infrared oscillations is about five times 
greater than that observed at radio wavelengths. The mean 
radio - infrared spectrum for the four-day period illustrated 
in Fig 3 is shown in Fig 4. 

A striking feature of the light curves in Figs 1 & 2 is 
the infrared oscillation which starts around MJD 51318.59 
which does not appear to have a mm counterpart, unlike the 
other four simultaneously observed events. We have care- 
fully checked the data reduction techniques to see if this 
was due to human error, but found no evidence of this. We 
have no clear physical intepretation of this phenomena, ex- 
cept to suggest that it was due to strong and variable ab- 
sorption which only signficantly affected the lower- frequency 
emission (for example the optical depth to free-free absorp- 
tion, r cx v~ 2,1 ). As noted above, the overall radio-mm- 
infrared spectrum at the time of these observations (Fig 4) 
was steeper than previously observed, perhaps also indica- 
tive of some absorption. If this is a correct interpretation of 
the 'failed' mm event, then we would have expected the ra- 
dio emission to have been completely absorbed at this time 
also. 



4 ENERGETICS AND MASS OUTFLOW RATE 

The radiative luminosity of these oscillations is large - for a 
fiat spectrum of amplitude 300 mjy from 1 GHz to 1.4 x 10 5 
GHz (= 2.2/xm), at a distance of 11 kpc, it is 3 x 10 37 erg s _1 
(the time-averaged radiative luminosity is around half this 



value) . As is the case for all synchrotron emitting plasmas for 
which adiabatic expansion losses dominate, this is likely to 
be a significant underestimate of the power being supplied 
to the jet. Furthermore we assume the emission arises in 
a partially self-absorbed jet which retains the same power- 
law distribution of electrons (i.e. p = 2.6 where N(E)dE oc 
E~ p dE) as observed in optically thin ejections (F99). In 
this situation the flat spectrum is produced by a conical, 
partially self-absorbed jet (e.g. Blandford & Konigl 1979; 
Reynolds 1982). 

In addition several factors which can further affect the 
energy budget are uncertain, in particular whether or not 
the small ejections share the same bulk relativistic motions 
as the larger ejections, whether each radiating electron has 
an associated cold proton, and what the filling factor (ie. 
the effective volume) of the ejecta is. The procedure for cal- 
culating the energy and mass of the ejections is as follows: 

• Transform observed flux densities and frequencies back 
to their rest frame (identical if no bulk relativistic motion). 

• Integrate rest-frame luminosity. 

• Calculate maximum emitting volume - in this case 
based on the five-minute rise time this is 3 x 10 39 cm 3 . The 
effective emitting volume is this volume multiplied by a 'fill- 
ing factor', /. 

• From the volume, spectrum and luminosity, calculate 
cquipartition magnetic field, and corresponding minimum 
internal energy. 

• For baryonic case, add one proton for each electron. 

• For cases with bulk relativistic motion, add in kinetic 
energy and multiply by two, under the assumption that ob- 
served emission was dominated by one (approaching) com- 
ponent only. 

• Divide by repetition quasi-period of oscillations to ob- 
tain time averaged energy and mass outflow rate. 

We have tabulated results for different cases in table 1; 
for bulk relativistic motion we have used the Doppler factors 
corresponding to (3 = 0.98, 6 = 66° from F99. 

A significant constraint is that the lack of evidence for 
synchrotron losses (based on the similarity of decay rates 
at widely different wavelengths) at 2.2/xm on a timescale of 
~ 10 min, implies that B max 30 G. For bulk motions with 
Doppler factor S, this limit is shifted slightly to Bmax^ -1 ^ 3 
which, for S = 0.34 in this case means the limiting field is 
~ 40G, instead of ~ 30G, not a major difference. Thus while 
reducing the emitting volume via the filling factor (see Ta- 
ble 1) decreases the minimum energy, the stronger derived 
cquipartition magnetic field is irreconcilable with the ob- 
served minimum lifetimes. As already noted in F97, a field 
of order 10 G will cause a cut-off in the spectrum of the os- 
cillations around the optical band. Reducing the magnetic 
field below the equipartition value results in the internal 
energy being dominated by the electrons, for which total 
energy E oc B~ 3 / 2 . Because of this constraint, the realistic 
minimum energy cannot be reduced much below the value 
for a non-relativistic non-baryonic ejection with filling fac- 
tor / = 1.0, which is 4 x 10 41 erg, with a corresponding 
time-averaged power requirement of 3 x 10 38 erg s _1 . This 
corresponds to a radiative efficiency for the outflow of < 5%. 
This minimum power requirement is not strongly affected by 
our assumption of the spectral form of the electron distri- 
bution (e.g. for p — 2.0 the minimum power is reduced by 
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Figure 3. Ryle Telescope observations of GRS 1915+105 before and after the JCMT/UKIRT observations. Each of the top panels 
corresponds to 0.2 days, and their distribution with respect to the JCMT/UKIRT observations is indicated in the lower panel. It is clear 
that strong radio oscillations were occurring before and after the JCMT/UKIRT observations, with a similar quasi-period, but lower 
amplitude. In fact the RT observations show that the oscillations were probably continuous between MJD 51309-51320. 
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Table 1. Calculation of radiative luminosity, cquipartition magnetic field, total energy and jet power & mass-flow rate for the oscillations 
reported here, given different physical assumptions. T is the bulk motion Lorcntz factor, / is the 'filling factor'. In these calculations a 
distance of 11 kpc and Doppler factors for relativistic bulk motion which are the same as those reported in F99 are all assumed. Mass 
flow rate Mj e t an d jet power P are based upon one ejection every 20 min. For more details, see main text. 



only a factor of three). Only in the case of baryonic ejec- 
tions at high velocities does the large number of lower en- 
ergy electrons become significant, as each has an associated 
proton. As a result it does not matter whether the emis- 
sion arises in a discrete 'plasmon' or an internal shock in a 
quasi-steady flow. It is interesting to compare the power for 
baryonic ejections with T — 5 , 2 x 10 43 erg s _1 , with that 
calculated for the same criteria for the 'major' ejections in 
F99, 2 x 10 39 erg s~\ This is due to the observed optically 
thin cm spectrum, S v oc v~°' s , being assumed in F99 to 
have no high-frequency excess, and therefore a much lower 
integrated luminosity than the flat spectrum oscillations re- 
ported here. Further prompt mm and infrared observations 
during 'major' outbursts are required to investigate this. 



5 DISCUSSION 

M98 have shown that the wavelength-dependent time de- 
lays (radio-radio and radio-infrared) observed from GRS 



1915+105 can be approximated by a 'van der Laan' (1996) 
model for an expanding plasmon. However, as discussed in 
FP98 such a model does not well describe the observed flat 
spectrum, which seems instead to be better modelled by a 
partially self-absorbed conical jet of the type developed for 
AGN (e.g. Blandford & Konigl 1979; Reynolds 1982). 

More recently Kaiser, Sunyaev & Spruit (2000) have fur- 
ther applied a internal shock model to the major 1994 radio 
outburst of the source reported in MR94. In their model 
they require approximately the same amount of energy to 
be associated with the events as calculated for a plasmon 
model in MR94, but the power requirement is less as they 
spread the energy input over a much longer period. However, 
with repeated oscillations as observed here, this cannot be 
the case, where an entire accretion - ejection cycle is re- 
peated on the timescale of ~ 20 min which we have used to 
calculate P and Mj ct in table 1. Therefore this model cannot 
be used to evade the enormous amount of continuous power 
required to generate the observed repeated ejection events 
(this is not an argument against their model, but one against 
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Figure 4. The mean radio-mm— infrared spectrum of GRS 
1915+105 in the interval MJD 51316-51320, centred on our si- 
multaneous infrared/mm observations. The spectrum appears to 
be significantly inverted at longer wavelengths. Vertical bars on 
each point are the statistical standard deviation for the data set, 
and reflect the relative amplitude of observed variability at each 
wavelength (not measurement errors). 



using it to evade the huge power requirements). Importantly, 
unless (a) there is a bright mm-infrared contribution from 
the large ejections which has not to date been observed, and 
(b) it is only the spatially resolved ejections (F99; RM99) 
which have bulk relativistic motion and a baryonic content, 
then GRS 1915+105 injects more energy and matter into 
the outflow during periods of repeated small events than it 
does during the large ejections. 

Note also that Belloni, Migliari & Fender (2000) have 
found that jet power, calculated as above, appears to be 
anticorrelated with accretion rate as inferred from X-ray 
spectral fits, for a small sample of observations with quasi- 
simultaneous infrared and X-ray coverage. 



6 CONCLUSIONS 

We have reported giant repeated oscillation events from the 
black hole system GRS 1915+105 observed simultaneously 
at mm and infrared wavelengths. Contemporaneous radio 
observations indicate that these observations were near the 
end of a sequence of i> 10 days of oscillations with ~ 20-min 
quasi-periods (i.e. <; 700 discrete ejection events). We have 
investigated in depth the energy and mass flow associated 
with such events, seeking to minimize the very large required 
power. However the magnetic field has an upper limit im- 
posed by the lack of observed radiation losses in the infrared 
band, so that reducing the effective volume by means of a 
small filling factor cannot significantly reduce the required 
power. Given the repeated nature of the events, and their 
almost certain association to an accretion cycle, we cannot 
spread the minimum required energy over a longer timescale 
than the repetition quasi-period. As a result we find that at 
least 3 x 10 38 erg s _1 is required to be channelled into the 
formation of the cjecta; a very significant fraction of the ac- 
cretion energy unless the black hole is very massive indeed. 
The enormous energy budget associated with baryonic ejec- 



tions at bulk relativistic velocities may rule out this possi- 
bility, although at present we do not know which is more 
unlikely: no baryons or low velocity ejections. 
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